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Abstract 

We propose a stochastic map model of economic dynamics. In the last decade, 
an array of observations in economics has been investigated in the econophysics 
literature, a major example being the universal features of inequality in terms of 
income and wealth. Another area of inquiry is the formation of opinion in a society. 
The proposed model attempts to produce positively skewed distributions and the 
power law distributions as has been observed in the real data of income and wealth. 
Also, it shows a non-trivial phase transition in the opinion of a society (opinion 
formation). A number of physical models also generates similar results. In particular, 
the kinetic exchange models have been especially successful in this regard. Therefore, 
we compare the results obtained from these two approaches and discuss a number 
of new features and drawbacks of this model. 



1 Introduction 



It is known that the distributions of income and wealth possess some robust 
and stable features which are independent of the economy-specific conditions 
\T\. But the exact form of the distribution is still debated [2]. It has been a 
tradition in the economics literature to model the left tail and the mode of the 
distributions of the incomes with a log-normal |3] distribution and the right 
tail with a Pareto distribution i.e., a power law [1]. However, a number of 
recent studies in econophysics literature argue that the left tail and the mode 
of the distribution is best described by the exponential or gamma distribution 
and the right tail of the distribution follows a power law (See Ref. [5|6|8] ). Also, 
a significant number of attempts have been devoted to explain the emergence 
of a consensus in a society and the possibility of a non-trivial phase-transition 
in the average opinion of the society (See Ref. 
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Here, we propose a stochastic map which attempts to produce all of the afore- 
mentioned features at different limits and in addition, it shows some new 
features as well. Below, we propose the model and discuss its relative merits 
and demerits with respect to the kinetic exchange models which also generates 
similar results. For related literature, see Ref. pjr|l8|IT5] for theoretical and nu- 
merical results on the kinetic exchange models of markets. See Ref. [9|10llll] 
for detailed analysis of the limit distributions of the solution of random dif- 
ference equations (of the form x{t) = a{t)x{t — 1) -|- we will use only a 
few very particular instances of it). Ref. [12] characterizes gamma distribution 
which arises from a random difference equation. Ref. |2DJ was the first paper 
that suggested that multi-agent interactions in the kinetic wealth exchange 
models can be simplified and and viewed from from the point of view of a sin- 
gle agent. Ref. [T^ was possibly the first paper that connected the stochastic 
maps and the kinetic exchange models. 



2 The map 

The Kolkata kinetic wealth exchange models have been successful to generate 
a realistic description of the income/wealth distributions (See ref. [S]). Many- 
body interactions are the key ingredients of this class of (kinetic) exchange 
models. In this paper, we mainly use a single agent framework to discuss the 
distributional issues. There is no agent-agent interaction in this model. In this 
regard, our treatment is closer to the representative agent paradigm of the 
modern economics. The specific structure of the market is also not discussed 
here. We represent the interaction of the agent with the market with a map 
of the following form. 



where Ai, A2 and A3 are linear functions of a single parameter A with < 
A < 1 such that < Aj < 1 for i=l, 2 and 3. We assume that —00 < 
n < 1, 9 takes value either a positive, finite value (we assume 9 = 1, for 
convinience) or 00 (or a sufficiently large value) and e^, ~ uniform[0, 1] and 
independent, unless specified. In all the simulations, we have assumed m(0) 
= 1. We denote the time index by subscripts for exogenous random variables 
(e and ^) which we shall drop when no confusion arises. While the parameter 
space is seemingly too large to be considered in details, we shall restrict it 
considerably by assuming very simple forms of Aj for all i. 

We can interpret it in the following way. An agent has m(t) amount of wealth 
(or money) at time t of which he saves a random fraction and invests the rest. 
His return from investment is represented by the additive term ^tAg. However, 



m(t) 



min {(Ai + etX2)m{t - 1) + ^A^, 9} 
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there is an upper limit of m{t) represented by 6. Since we interpret the random 
multiphcative term as the savings propensity, we assume that Ai + €1X2 < 1 
for all t i.e., savings propensity is never greater than 1. Our interest lies in 
finding the pdf (probability density function P{m)) in the steady state. 

It may be noted that the above map (ignoring 9) has the general form 

m{t) = a{t)m{t - 1) + h{t) (2) 

which has been studied in great details by Ref. pilUfllj . In this paper, we have 
focused on a few particular instances of it for our purpose. Ref. [2] mapped 
the asset exchange models into 'random iterated function systems'. But it 
was concerned with the 'yard-sale' model and the 'theft-and-fraud' model (see 
Ref. [13]) whereas we focus on the CC and CCM models (the models that 
introduced fixed and distributed savings propensities in the kinetic exchange 
models; see Ref. [5]) and that is why we have borrowed the basic structure (a 
fraction of the total wealth is saved and the rest is invested) from the Kolkata 
wealth exchange models (i.e., the CC and CCM models). 



3 Statistical features of the map at some limits 

3.1 A positively skewed distribution 



We assume Ai = A, A2 = A3 = 1 — A, ri, = 1 and 6 — t- 00. The equation becomes 
m{t + 1) = {A + e(l - A)}m(t) + ^(1 - A). (3) 

See Fig. [T] for steady state distributions corresponding to different values of 
A. We restrict our attention to the case n = 1 as opposed to other positive 
numbers, for two reasons. One, this choice of n makes the average a constant, 
independent of n (as the CC and CCM models; see Ref. |8]) and secondly, to 
simplify the calculations. We will relax this assumption in Section 13.21 where 
we will assume n is negative and we will see that the exponents of the power 
law distributions differ for different values of n. 

Though we do not know the exact algebraic form of the steady state distri- 
butions produced by Eqn. [3l we can find out its moments in order to describe 
them qualitatively. We can ignore the time index in the steady state. Taking 
expectations over both sides of Eqn. [21 we get 

(m) = 1. (4) 
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Fig. 1. 

Steady state distributions generated by Eqn. [S] Three cases are shown above 
viz., A = (+), A = 0.4 (x), A = 0.7 ( *). All simulations are done for ~ 10' 
time steps and averaged over for ~ 10^. 

Also, we have the variance of m as 



Vim) = (x") - (x)' 



(5) 



where x = (A + e(l — A)) m + ^(1 — A). Now, we make a few almost trivial 
observations. Note the fact that (x) = 1. Hence, (m^) can be written as 
V{m) + (m)^ (by defn. of V{m)) i.e., V{m) + 1. Also, e and ^ are uniformly 
distributed. Therefore, (e) = 1/2 = and (e^) = 1/3 = (^2) (recall that V{e) 
= 1/12 = V{C,)). Using all of these and by expanding Eqn. |5]we get 



V{m) = [A^ + A(l - A) + i(l - A)2] (\/(m) + 1) 
+i(l-Af + i(l-A2)-l. 



Simplifying the above expression we get the result for A 7^ 1, 
V{m) 1^1-^^ 



2 \2 + A 



(6) 



which clearly shows that the distribution tends to a delta function as A — > 1. 
With A = 0, Eqn. [3] produces a distribution with a sharp peak at m very close 
to 1. However, for A > 0.3 this sharpness goes away. 



One can consider an even simpler case with =^t for all t so that the map 



4 



becomes 



m(t + 1) = Am(t) + e(l - A) (m(t) + 1) (7) 



which is almost the same as the usual kinetic exchange model with a constant 
savings factor i.e., the CC model (See Ref. [2Z]; See also Ref. [S]), 



m,(t + 1) = \mi{t) + e(l - A) (m,(t) + mj{t)) . 



The variance of the distribution generated by EqnJT] is given by 

V(^n) ^ (i^) . (9) 



For the sake of completeness, we mention that the variance of the distribution 
generated by Eqn. [HI is given by (See Ref. [28] ) 

V(m) = ~ . (10) 
^ ^ (1 + 2A) ^ ' 



It is noteworthy that for A = 0, the distribution generated by Eqn. [7] has 
a peculiar form which is known as Dickman distiburtion (see Ref. [T^). For 
m < 1, P{m) is fiat and for m > 1, the distribution has a downward slope 
(see Fig. [2]). See Ref. [I6], for this type of maps which has been used in number 
theory, in biology (see Ref. [17]) and in many other areas (see Ref. [IS]). See 
Ref. [H] for simulations on the pdf generated by Eqn. [7] with A = 0. 

It is to be noted that we cannot generate an exponential distribution from 
Eqn. [7] since its maximum variance is 1/2 (for A = 0) whereas the variance of 
the exponential distribution generated by the following Eqn. (See Ref. [8]) 

mi{t + l) = e{mi{t) + mj{t)) (11) 



is unity. Hence this model can not generate purely exponential distribution 
which is a drawback since recent studies argue that the income/wealth distri- 
butions in the real world fits very well with exponential pdfs [T|5|6] . It is trivial 
to note that in the usual kinetic exchange models, it is the presence of the 
trading partner's wealth mj{t) (in the i-th agent's wealth evolution equation; 
Eqn. [TTj) that contributes to the higher variance. 
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Fig. 2. 

Steady state distributions generated by Eqn. [7] are plotted in the semi-log 
plot : Three cases are shown above, viz., A = (+), A = 0.4 (x), A = 0.7 
( *). All simulations are done for ~ 10^ time steps and averaged over for 
~ 10^. Inset : The same data set in the usual plot. 

3.2 Power law distribution 



We assume Ai = A, A2 = A3 = 1 — A, —00 < n < and 6 ^ 00. 
The relevant equation is 

m{t + 1) = {A + e(l - A)}m(t) + ^(1 - A)". (12) 

Ref. [22)123] studied this type of discrete stochastic equations (see Eqn. [2]) as a 
generic model for generating power law pdf. But a necessary condition for the 
mechanism to generate a power law pdf is that the random multiplicative term 
a{t) (A + e(l — A) in Eqn. [T2|) must be greater than 1 sometimes. In the current 
context, this condition is not satisfied (according to our interpretation, (A + 
e(l — A)) is the savings propensity and hence it is always less than unity). To 
avoid this problem, we assume a population of agents each of which interacts 
with the market according to Eqn. [T^l Note that the agents do not interact 
among themselves. In the spirit of the kinetic exchange models, we show that 
if there is a population of N agents with different A but each of the income 
evolution process is modelled by Eqn. [121 then a power law in the income 
distribution will be observed. However, there is an important difference with 
the usual kinetic exchange models which are completely conservative. Any 
trading activity in such markets would be a zero — sum game i.e., if somebody 
wins then his/her trading partner has to lose. We relax that assumption here. 
For simulations, we assume that there are 200 agents each endowed with an 
initial wealth equals to 1. Each of the agent's wealth evolution is governed by 
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Fig. 3. 

A power law distribution {n = —20 in Eqn. [T2|) . The straight hne with slope 
-1 is drawn as a guide. However, here we take 0< A <0.33 because otherwise 
the average value of m becomes too large; e.g., Eqn. [13] shows that for A = 

0.5, (m) = 2^1. 
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Fig. 4. 

A power law distribution (n = in Eqn. [T2l) . The straight line with slope -2 

is drawn as a guide. 

Eqn. [121 Clearly, there is no interactions between the agents. However, the 
agents are assigned different As which are fixed over time for any given agent 
(see Ref. |29]). In particular, we assume that A is uniformly distributed among 
the agents. For each agent, simulations are done for ~ 10^ time-steps and the 
corresponding pdfs are averaged for ~ 10^ time-steps. Resulting distributions 
are averaged over all agents. 

Following Ref. [20] , a very simple proof is considered below. Note that (in the 
steady state) by taking expectations on both sides of Eqn. [T21 we can rewrite 
it as 

(l-A)^-"(m) = l (13) 
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Fig. 5. 

The pdf of wealth distributions are drawn for the case n = —20 with global 
savings propensities 0, 0.1, 0.2, 0.25, 0.3 (from the left to the right) in the 
log-log plot. Also, we've drawn in the same diagram (the dotted line). 




Fig. 6. 

The pdf of wealth distributions are drawn for the case n = with global 
savings propensities 0.1, 0.5, 0.8, 0.9, 0.95, 0.98 (from the left to the right). 
Clearly, the average values of wealth at any given savings propensity are very 
close to the theoretical values predicted by Eqn. [131 Also, we've drawn m~^'^ 
in the same diagram (the dotted line). Inset: The same diagram in log- log 

plot. 



By taking total differentiation and rearrenging terms, we get 



dm 



-m 



n 



(14) 



where m represents (m). Hence, the average amount of money held by an 
agent with a particular A is given by Eqn. [T31 Also, the relation between the 
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Fig. 7. 

The pdf of wealth distributions are drawn for the case n = with maximum 
savings propensities Xm = 0.9, 0.95, 0.975, 0.99 (from the left to the right) in 

the log- log plot. The savings propensities of the agents are chosen 
deterministically following the rule that Aj = {i/N)XM- Also, is drawn 

for reference (the dotted line). 

distribution of A (i.e., p(A)) with that of m is given by the following Eqn. 
P{m)dm = p{X)d\. (15) 




Eqn. [13] and [15] shows that in an population of agents with uniformly dis- 
tributed A, the distribution of m would be 

n-2 

P{m) ra "-i . (16) 



It is trivial to note that the CCM model (See Ref. |29]; See also Ref. [8]) 
derived the result analogous to the case where n = (see Fig. Hj). Also, if n is 
large in modulus then the steady state would be power law distribution with 
power -1 (see Fig. [3]) (we take n = —20 just for expository purpose). 

Following Ref. [S], we can support the above result for the case n = 0. Let 
us rewrite Eqn. [T2] with the agent index as the following 

m,{t + 1) = {\i + e,(l - \i)}m,{t) + ^^(1 - A,)"- (17) 

which can be rewritten (using Eqn. [T3il as 

rriiit + 1) = {Xi + q(1 - Xi)}mi{t) + ^^(1 - Xi){mi). (18) 

The above equation holds true for all agents i. Since we focus on average 
money holding at different savings propensities, we randomly pick any two 



9 



agents i and j (with savings propensities Aj and \j respectively) and multiply 
their corresponding wealth evolution equations (fTSjl . to get 

mi{t + l)mj{t + 1) = fljmi{t)mj{t) + fij{mi{t))mj{t) 

+4m,(t)(m,(t)) + 4(m,(t))(m,(t)) (19) 

where f^s are functions of e, ^, Aj and Xj for all /c. We approximate each of 
the quadratic quantities by a mean quantity m^. Note, that here we assumed 
that rrii can be replaced by its average value {rrii), which holds true only if its 
variance is small and that requires n to be as small in modulus as possible in 
Eqn. [T71 We consider n = 0. Then, we have the following equation 

m'^it + 1) = a(t + l)m^(t) where a is a function of f'^s. . (20) 

Ref. [2T|8] studied this map and showed that as t — )■ oo, 

P(m) ~ m-^. (21) 



Note that this argument holds only if n = 0. 

We consider the cases where all agents have the same savings propensities 
which we call the global savings propensities (this case is almost identical 
to the one described in Section 13.11 in that we consider agents who have the 
same A; the difference is that there we considered only one agent whereas here 
we consider identical agents. Also, another obvious difference is that here 
n < 0). The pdfs are drawn for individual savings propensities in Fig. Oand 
El Let m have a distribution /(m) in the steady state (given A). With a slight 
abuse of notation, we write the value of /(.) evaluated at m = m as fifn). 
Then Fig. [5] shows that (for n = —20) given A at m = (m), f\{{m)) ^ (m)^^. 
Similarly, for n = given A at m = (m), fx{{m)) ^ (m)^^'^. This shows that in 
the case where n is large in modulus, the fx{m) curve is a good approximation 
of the average of all the individual wealth distributions which is not true in 
the other case where n = 0. The reason is that the probability of finding an 
agent within a specified interval of wealth decreases very rapidly with increase 
in wealth in the case where n = compared to the other case (this can be 
checked very easily from Eqn. [T^ . Hence, though the uppermost points form 
a locus of fx{{m)) ~ (m)""'^ for n = 0, the average distribution falls much 
more rapidly following P{m) ~ m~^. 

We have also studied the effects of the maximum A present in the population 
on the steady state distribution (for the case n = 0). Following Ref. [21], we 
the assign the savings propensities to the agents deterministically, according 
to the rule that the i-th agent's savings propensity Aj = {i/N)\M where N 
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is the number of agents in the population, \m < 1 is the maximum savings 
propensity in the population. See Fig.[7]for the results. Curves from left to right 
represent the distributions for increasing Xm- Clearly, as Am decreases, the 
power law interval shrinks. Therefore, higher savings propensities contribute 
to the formation of the power law interval. This conclusion is consistent with 
the findings in Ref. [2lj. We should also mention that according to Eqn. [T3| 
higher A implies higher average wealth at that particular A. This fact also 
supports the findings in Fig. [3 

3.3 A non-trivial phase transition 

Here, we briefly mention that a special case of the map shows a non-trivial 
phase transition. If we assume Ai = A2 = A, A3 = 0, n = 1 and = 1. The 
equation becomes 

m{t + l) =min{\{l + e)m{t),l} (22) 

This map has been used to model opinion formation. Assume that an agent can 
have an opinion m{t) within and 1 at any time t. A is the agent's conviction 
parameter. After each interaction with the society, his opinion changes by a 
random fraction. However, the maximum value of opinion that he can have, 
is 1. Ref. considered this map and showed that this map shows a phase 
transition (with respect to A) in the average value of m. The critical value 
of A found from simulations, is Ac ~ 0.68 and it only mentions an analytical 
result that Ac ~ 0.6796. For numerical details on this particular map, see Ref. 
[26] . Interestingly, this map has a parallel in the kinetic exchange models, also 
discusses in the same reference (which has a critical point at Ac ~ 0.6667). For 
detailed analysis on the nature of such transitions see Ref. [5U|^ . 

Let us now focus on the cases with A < 1. Consider the map 

m(t + 1) = A(l + e)m(t). (23) 

Let us ignore m{t) for alH < r where r is sufficiently large. For any time step 
T>r, 

m{T) = X^'^l + eT)(l + eT-i)...(l + eT-.+i)m(r). (24) 

Following Ref. ^26j, we can argue that for A < Ac, the product of all the terms 
multiplied to m(r) in the above equation is less than unity on average. But 
this term tends to 1 as A — )■ Ac. Therefore at the critical point (A = Ac), by 
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taking logs on both sides we can write 

-log A, = -^f]log (1 + e,). (25) 

^ k 

As T — )■ oo, we apply the LLN and the r.h.s. of the above Eqn. converges to 
(log (1 + e)). By applying Jensen's inequality, we see that 

-log Ae<log((l + e)) (26) 

i.e., Ac > 2/3. Therefore the critical point of the map model is greater than 
that of its kinetic exchange version. By numerical calculations from Eqn. [251 
we see that Ac ~ 0.67954 which is very close to the value reported in Ref. [26]. 
One could also find that (log (1 + e)) = 2 log 2 — 1 and therefrom find Ac using 
Eqn. [23 (see |26|). 

It is trivial to note that for A < Ac, m(t) tends to 0. Also at A = Ac, the 
distribution has an abrupt change in variance as has been shown numerically 
in Ref. [26]. Here, we present a very short proof. By Eqn.[24l (m(T)) = (m(r)) 
at A = Ac and hence for T — r sufficiently large, 

m(r) - (m(r)) ^ Aj--(1 + eT)(l + eT-i)...(l + eT-.+i)(m(r) - (m(r))). (27) 



By squaring and taking expectations on both sides, we get 



V{m{T)) ^ {X'f-^\l + erYi^ + ^T-if-il + er-.+i)')\/(m(r)) 
> (Aj-^(1 + eT)(l + er-i)...(l + er-.+i)) V(m(r)) 
= V{m{r)). (2J 



To get the last inequality, note that V{x) = (x^) — {x)^ > for all x (unless 
it is the trivial case that x = (x)). Therefore, V^(m(T)) > Z(^T~T)V{nT'{'^)) with 
Zk > 1 for all k and for all r. Hence, as t — ?■ oo, V{m{t)) — t- oo at A = Ac. 

To illustrate the problem of non-stationarity, we focus on another particular 
case where A — )■ 1, e ~ uniform[0, e] with e — > such that (e) = — log A i.e., 
(e) = — log (1 — (1 — A)) ^ 1 — A. From Eqn. [21|, by taking logs we get 

T 

log m(T) = (T - r) log A + ^ log (1 + e^) + log ni{T) 

k 

T 

~E(efc-(efc)) + log m(r). (29) 

k 
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Multiplying both sides by 1/yT and by applying CLT (assuming T — )■ oo), 
we get 

^ log m(T) A^(0, (T^) in distribution, (30) 



where o"^ is the variance of e. Hence, m{T) is distributed log-normally. But 
evidently, the distribution is not stationary. 



3.4 The law of proportionate effect 



We assume Ai = A2 = 1, A3 = and — ?■ 00. Note that this is the limit of 
the case discussed above with A = 1 and 6 sufficiently large. In this case, the 
equation becomes 

m(t + l) = {l + et+i)m{t) (31) 

which is very wellknown as a generator of log-normal distribution (See Ref. 
[3]). To find the distribution of m(t), note that the above equation can be 
rewritten (by iteration) as 

t 

m{t + l) = X{{l + ek)m{l) (32) 

k 

and by taking log on both sides, it can be written as 

t 

log m{t + 1) = X! (1 + ^k) + log m(l) 

k 

t 

oi'^ek + log m(l) (for e very small). 

k 

(33) 



By applying the Central Limit Theorem, we see that log m{t) is distributed 
normally, hence m{t) is distributed log-normally. However, evidently the dis- 
tribution generated this way is not stationary either. 
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4 Summary 



In recent years, a number of economic regularities have been investigated in 
the econophysics hterature, one of the most prominent themes being the distri- 
butions of income/wealth. The two candidate distributions for explaining the 
left tail and the mode of the distributions of income/wealth are Gamma and 
Log-normal. There is a consensus that the heavy right tail is best described 
by a power law (See Ref. [HH]). Also, the process of opinion formation due to 
interactions among numerous agents have been studied in details. The phe- 
nomena of emerging consensus and a phase transition in the opinion formation 
have been tried to be modelled in many ways (See Ref. [2H|26j ). 

In this paper, we examine a stochastic map to model the economic process of 
income/wealth distribution and the sociological process of opinion formation. 
We showed that this particular map can generate a positively skewed pdf and 
a power law pdf. It shows a non-trivial phase-transition and in another limit, 
it coincides with a very well known generator of log-normal pdf. However, 
there are several drawbacks of this approach. One is that the algebraic form 
of the pdf generated by the Eqn. [3] is not known and it is not a F pdf (at 
one limit it has a sharp peak which eventually goes away). Second is the well 
known problem with Eqn. [31] that the distribution generated by this equation 
is not stationary. Third, purely exponential distribution cannot be generated 
by this model. But there are a number of recent studies concluding that the 
income/wealth distribution can be modelled by exponential pdf p!|5|6] . Pre- 
viously, the kinetic exchange models have been successful to model many of 
these phenomena pi26] . We show that there are some new results produced 
by this map model whereas other results mostly conform with those derived 
in the kinetic exchange models. 
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